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DYNAMICS OF SCREW DISLOCATIONS: A GENERALISED 
MINIMISING-MOVEMENTS SCHEME APPROACH 

GIOVANNI A. BONASCHI, PATRICK VAN MEURS, AND MARCO MORANDOTTI 


Abstract. The gradient flow structure of the model introduced in |CC99) for the dynamics of screw 
dislocations is investigated by means of a generalised minimising-movements scheme approach. The 
assumption of a finite number of available glide directions, together with the “maximal dissipation 
criterion” that governs the equations of motion, results into solving a differential inclusion rather than 
an ODE. This paper addresses how the model in |CG99| is connected to a time-discrete evolution 
scheme which explicitly confines dislocations to move each time step along a single glide direction. 
It is proved that the time-continuous model in |CG99| is the limit of these time-discrete minimising- 
movement schemes when the time step converges to 0. The study presented here is a first step 
towards a generalization of the setting in |ACS081 Chap. 2 and 3] that allows for dissipations which 
cannot be described by a metric. 


Keywords: Motion of dislocations, generalised gradient flows, minimising-movement scheme, energy dissipa¬ 
tion inequality. 
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1. Introduction 

This paper is devoted to the application of the minimising-movement framework to the model 
introduced in |CG99) for the evolution of a finite number of screw dislocations, under the constraint 
that dislocation movement only occurs along a finite set of glide directions. This constraint enforces 
that defects only move along the so-called slip planes, which are determined by the crystallographic 
structure. The direction followed by each dislocations is the one dictated by a maximal dissipation 
criterion: it is the glide direction which is closest in direction to the Peach-Kohler force that acts on 
the dislocation. 

Yet, this confinement of the dislocation motion along glide directions is not captured by the model 
in |CG99| . Indeed, this model allows for dislocation glide along more complex paths: dislocations 
can switch from one direction to another one, or even move along curved lines. These behaviours are 
called cross-slip and fine cross-slip, respectively, in |CG99) . Roughly speaking, at the turning point 
in cross-slip the dislocation switches direction from the most dissipative one to another one which 
becomes equally dissipative. This switch between one glide direction and another happens at a much 
faster time scale in fine cross-slip, giving the dislocation a curved trajectory. 

The aim in this paper is to understand how such curved trajectories emerge from dislocation 
motion along glide directions. Our main question is: 

How does the evolution model in |GG99) connect to a different model which confines 
dislocations to move along glide directions only? 

Our main contribution is to solve this question by proposing a time-discrete variational model (i.e. a 
minimising-movement scheme) for the movement of dislocations, which confines the motion of any 
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dislocation at each time step to a single glide direction, and by proving that the limiting equation, 
as the time step goes to zero, is given by the model in |CG99| . Interestingly, the conhnement of 
dislocations to move along glide directions leads to a dissipation which cannot be described as the 
square of a distance (it will not satisfy the triangle inequality). Therefore, the well-known convergence 
of minimising-movement schemes (see e.g. [AGSOS) ! does not apply, and our result is a hrst step in 
generalizing the convergence of minimising-movement schemes. 

After a brief prologue on dislocations in Section |1.1[ we introduce the model from |CG99) in 
Section [I| an d our modihed setting in Section HI We then provide an intuitive example (see 
Example 1.1) to illustrate the dynamics of screw dislocations. In Section 1.5 we describe and discuss 
Theorem 3.1 which is our main result that describes the connection between our modihed setting 


and |GG99 


1.1. Screw dislocations. Dislocations are line-defects on the atomic length-scale in the crystallo¬ 
graphic structure of the material. Typical metals contain many dislocations (for example, cold-rolled 
metal has a dislocation density of 10^^ m^ |HBni| p. 20], which translates into 1000 km of dislocation 
line per cubic millimetre). Their collective motion results in plastic deformation of the material on 
length-scales between 1 /rm and 1 mm. Figure shows two straight segments of a dislocation line 
in a cubic crystallographic lattice. The result of dislocations being line-defects is that they induce 
a stress in the material. Dislocations elsewhere in the material may move as a consequence of the 
stress. This results in an intriguing and complex system of interaction line-defects. For more details 
we refer to IhMIhboT] . 



Figure 1. Pictorial representation of edge (above) and screw (below) dislocations. 

(Source: https://commons.wikimedia.org/wiki/File:Vector_de_Burgers.PNG. 

Image by David Gabriel Garcia Andrade.) 

We study the evolution of a set of n straight and parallel dislocation lines of ’screw’-type (see 
Figure]^ in a three-dimensional elastic material undergoing antiplane shear. The body is modelled 
by an inhnite cylinder D x M. Therefore, we characterise the positions of the screw-dislocations by 
points in the cross-section D C Following the model proposed in |GG99j . we assume that the n 
screw dislocations zi,..., z„ G D are constrained to move along glide directions. The glide directions 
are determined by those directions along which the atoms in the material are most densely packed. 
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For instance, for a cubic atomic lattice, there are four glide directions given by {±ei,±e2} C (ej 
being the standard basis vectors in M^), and for a body-centred cubic or face-centred cubic lattice, 
there are six glide directions which span the triangular grid of equilateral triangles. We consider a 
more general setting in which the set of G 2N glide directions is given by 


G ■= {gi,-- 

•,gjv} c 

(1.1) 

which satisfies the basic properties 

g€G => -g&G, 

and span^ = 

(1.2) 


Under the maximal dissipation criterion described above, the velocity vector of any smooth solution 
t I—7- Z(t) = (zi,... ,Zn){t) is aligned along one of the glide directions at each time t. It turns out 
that generically, solutions of this type may not exist; this phenomenon is illustrated in Example o 
below, and the simulations in |BFLM1^ Sec. 4] confirm this observation. In order to obtain existence 
of solutions, one is forced to allow the motion of dislocations along a direction different than any of 
the glide directions, called fine cross-slip |CG99) (or ‘sliding motion’ as in the theory of discontinuous 
differential equations (see e.g. |dBBCK08| lb Note that fine cross-slip seemingly contradicts the model 
assumption that screw dislocations only move along glide directions. The variational evolution model 
that we propose here shows how fine cross-slip can be incorporated in the equations of motion, as 
the limit of rapidly alternating directions. 


1.2. The model in |CG99| . Before defining the evolution model which confines dislocations to 
move along glide directions only, we introduce the time-continuous model in |CG99) in more detail. 
Based on |CG99[ (1.10)], in |BFLM15] the following differential inclusion is posed to model the 
time-continuous dynamics of screw dislocations: 


dzi 

dt 


(t) G Fi{Z{t)), for all t > 0, i = 1,..., n, 


(1.3) 


where Z = (zi,..., z^)^ G O"' denotes the positions of n screw dislocations in a prescribed Lipschitz 
domain U C M^. Each dislocation is associated with a Burgers vector bj, which describes the direction 
of lattice mismatch due to the presence of the dislocation itself (see Figure]^. Due to the assumption 
of antiplane shear, in our case the Burgers vectors are oriented along the ea-direction, so that they can 
be described by a scalar quantity bi, which we call the Burgers modulus: bj = bie^. The quantities 
bi can be positive (5* = 1) or negative (6* = — 1) and their sign is responsible for the attraction or 


repulsion of dislocations, as it can be seen in (1.4) below. 


The set-valued function F)(Z) projects the force /j(Z) on the Uth dislocation onto the nearest glide 
direction. More precisely, we define the (nonlinear) multi-valued projection operator Pg : 'S? ^ 

by 


Pg^ ■= {(g-Og 


g G arg max g • O • 
gee J 


Then, FfiZ) = co Pgfi{Z), where co is the convex hull. 
|BM14I (4.4)], in which S takes the form 


The force /i(Z) = — Vzj£^(Z) is defined in 


S{Z) := ip{Z) -^'^'^-bibjlog\zi - zfi, (1.4) 

i=l j=l 

where ip G C'°°(D”') (see |BM14[ Lemma 5.1]) is bounded from above and satisfies (^(Z) —)• —oo as 
dist(zj,clD) —)• 0 for any i G N. The logarithmic interaction potential corresponds to the Peach- 
Kohler force induced by screw dislocations, and y? describes the effect of the traction-free boundary 
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condition at dQ. For our purposes it is enough to have S G \ S), where the set S of singular 

points is given by 

Z G S := dn^ U {Z G \ 3i ^ j:zi = zj}. (1.5) 

In |BFLM15) . local-in-time existence and uniqueness of solutions to ( |1.3[ ) is proven for suitable 
initial conditions. Here, a solution to (1.3) is defined to be an absolutely continuous curve which 


satisfies (1.3) for a.e. t > 0. The proof of well-posedness in |BFLM15) relies on the theory for 
differential inclusions developed by Filippov |FA88| . 

Next we discuss the result in |ADLGP15) in the context of our main question above. In |ADLGP15) 
a fully atomistic model is considered to describe the energy for a given configuration Z as a function 
of the atom spacing e. The dynamics are defined by imposing a dissipation potential which is the 
square of a norm which is minimal in the glide directions. The main results concern the derivation 
of the effective energy in the limit e —>■ 0 (which is also done in |ADLGP14) ). and the passage to 
the limit in the related minimising-movement scheme as the time step converges to 0. When the 
dissipation is chosen to be the square of the crystalline norm given by 

||x|| := inf |a -|- /3 I there exist a, /3 > 0 and g, g G ^ such that x = ag -|- /3g}, (1-6) 

then the evolution is obtained as a generalised gradient flow (see, e.g., [M ie m or Definition 

2.1). While this result connects rigorously a detailed atomic description of the dynamics of screw 


dislocations to the time-continuous evolution (1.3), it does not answer our question above. Indeed 


the dissipation in the atomic model is assumed to be derived from a norm, by which dislocations are 
not confined to move along glide directions. 

1.3. Evolution constrained along Q. To confine screw dislocations to move along glide directions 
only, we impose a discrete-in-time evolution model as a minimising-movement scheme. The special 
feature of this scheme is that the related distance is replaced by a ‘quasi-distance’ D given by 


d : X 


[0,oo], 


d(x,y) := 


|x-y|, ifx-yGM^, 
oo, otherwise, 


D : 


D2\n 


D2\n 


[0,oo], 


d2(x,Y) 


(1.7) 

( 1 . 8 ) 


2 = 1 


Since d is only finite along glide directions, D violates the triangle inequality. 


The minimising-movement scheme does not make sense for the energy £ defined in (1.4), because 


£ is not bounded from below on the set S of singular points (1.5). For this reason, we define a 
regularization of £ as follows: for any e > 0, we take G C'°°((M^)"') bounded from below, such that 
Tg = T on the closed set 

Ue := {Z G D”| dist(Z,5) > e}. 

To define the minimising-movement scheme, we take r > 0 as the time step, T > 0 as the end 
time, and Z® G D”' as the initial condition. A discrete-in-time solution Z^ G at the time points 
= kr is defined by 


rk+l 


G argmin<h(Z:^,X, r), /c = 0,..., [T/r] — 1, 

xer"' 


Z° = Z°, 

in which the functional <I> is given by 


(1.9) 




2r 























DYNAMICS OF SCREW DISLOCATIONS 


5 


By the definition of D, dislocations are confined to move along one glide direction only for each time 
step. 


1.4. Illustration of screw dislocation motion. Example o illustrates the evolution defined by 
(1.3). More involved examples concerning interacting dislocations are given in |BFLM15| . 


Example 1.1. This example is based on |CG99[ Figure 1]. We consider a single screw dislocation 
in a medium Q with glide directions given by ^ = {±ei,±e 2 } with e* the standard unit vectors. We 
impose a continuous force field E : —?• (0, oo)^. We choose F such that the set {x G | E(x) jj 

(ei + 02 )} equals the boundary of a bounded domain E, in which F ■ e .2 > F ■ ei, and outside of 
which F ■ ei > F ■ 02 . By (1.3), the dislocation will move along 62 if it is inside E, and along ei if it 
is outside M^. We show trajectories of the dislocation for several initial conditions in Figure]^ We 
limit our attention here to the direction of the movement and not to its speed. 



Figure 2. Typical setting of Example |1.H along with a few trajectories of the screw 
dislocation. The closed curve depicts the points x at which E(x) jj (ei + 62 ). 


The interesting part of the dynamics is the behaviour of the dislocation at dF. On dF the right- 


hand side of (1.3) is multi-valued. dF is called the ambiguity set. We distinguish the following three 


types of ambiguity sets: sources (dotted line), cross-slip (thin line), and fine cross-slip (thick line). 

If the screw dislocation is at a cross-slip point (say at the left-lower part of Figure]^ at t = ti, 
then it can move in any direction 8ei -|- (1 — 9)e2 for 6 G [0,1]. Whichever direction is chosen, at 
t = ti + At the dislocation is inside F for any At > 0 small enough, and hence it will move in 
direction 62 as soon as it enters F. 

Following a similar reasoning, at a fine cross-slip point (say at the left-upper corner of Figure]^, 
we conclude that a dislocation tends to move along a fine cross-slip set. If it hits the end point A, it 
will move along direction ei. 

Following a similar reasoning, we conclude that the initial-value problem may not have a unique 
solution at a source point. The word ‘source’ should be understood here in the sense that time paths 
will never cross it. 


1.5. Result and comments. Our main result is Theorem 3.1 It answers our main question by 


showing compatibility of the evolution constrained to glide directions with the evolution (1.3) that 
can cause dislocations to move in any direction. In particular, it gives a precise meaning to fine 
cross-slip in terms of oscillating choices for the preferred glide direction in consecutive time steps as 
the time step size converges to 0 . 
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More precisely, Theorem 3.1 states in what sense solutions to the minimising-movement scheme 


(1.9) converge to a solution of 


( 1 . 10 ) 


dz ■ 

-^{t) G coPg( - V , for all t E (0, T], z = 1,..., n, 

Zi{0) = z°, 

as the time step r —)• 0. By the definition of the regularization Tg, any solution of (1.10) satisfies 


(1.3) on the time interval [OjTg], where Tg is chosen such that any solutions Z(t) remains within 


the set rig. We notice that if we fix eq > 0, then Tgp > 0; moreover, it is easy to see that Tg is a 
non-increasing function of e. Therefore, e > 0 can be chosen arbitrarily. 


As a corollary of Theorem 3.1 we obtain existence of solutions to (1.3) up to the first time at which 


either two dislocations annihilate, or when a dislocation leaves the domain fi. This generalises the 


local-in-time well-posedness result of |BFLM15| . Theorem 3.1 also reveals how (1.3) can be written 


as a generalised gradient flow (see Section 2.3). 
|AnTXlP151 (3. 11)]. 


This gradient flow structure is also obtained in 


Theorem 3.1 is proved for all energies E E C'^((M^)"')nlT^’“((M^)"'). To the best of our knowledge. 


Theorem 3.1 provides the first extension of the theory in |AGSn8| Chapters 2 and 3] to dissipations 
which are not related to a distance. We wish to comment on further possible generalizations: 


Theorem 3.1 extends to higher dimensions for the particle positions with little modifications 
to the proof. We provide the details in Section [4.1| 

Our proof of Theorem 3.1 heavily relies on the regularity of the energy E. In view of the 
discontinuities in the projection operator Pg, it seems reasonable to require the force field 
—VE to be continuous for (1.3) to be well-posed. Hence, we do not aim to weaken the 
regularity conditions on E. 


By our regularisation Tg of T in (1.4) we cannot describe the annihilation of dislocations or 


the absorption of dislocations at dQ. In Section 4.2 we describe an open problem concerning 
the description of these effects in a variational framework. 

The remainder of this paper is structured as follows. In Section |2.3| w e define precisely the 
minimising-movement scheme (MMS) and describe the evolution model ^L3| in a variational frame¬ 


work by means of an energy dissipation inequality (EDI). Section 2.4 presents a priori estimates 


which are required for the proof of Theorem 3.1 which is contained in Section Section concerns 
the extension and the open problem related to Theorem 3.1 as mentioned above. 


2. Preliminaries 


2.1. Notation. Here we list symbols and abbreviations that we use throughout this paper: 


X 

Euclidean norm 


II ' II 

crystalline norm in with respect to Q 


• * 

dual norm of • 


x//y 

X is parallel to y, i.e. x/ x = y/ y 


CO 

convex hull 


d 

metric induced by 


d 

‘quasi-distance’, which is only finite along g E ^ 


b 

extension of the norm d to (M^)” 


D 

extension of the ‘quasi-distance’ d to (M^)"" 


s 

specific energy functional for edge dislocations; £ 


E 

generic energy functional; E E (^^((M^)”) Pi 

((M2)^) 


(^) 

(2.3) 


( 2 . 10 ) 

(T^ 

( 1 ^ ) 

(TW 

(1.4) 
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EDI energy-dissipation inequality 

<I> functional to be minimised in the D-MMS 

g, G glide direction in and set of all glide directions 

Ag cone around g 

MMS minimising-movement scheme 

Z particle positions; Z = (zi,..., G 

solution to the D-MMS at the fcth time step 
Zr step function related to (Z^)^; Z,- : [0,T] —)• (M^)” 
Z^ De Giorgi interpolant; Z^ : [0,T] —)• 


(gr 

(0I 

(^1 

Q 

(^1 

(2.81 


2.2. Glide directions. In this section we show several basic properties of the set of glide directions G 
defined in (1.1) satisfying the basic properties ( |1.2[ ). Figurej^illustrates an example. As a consequence 
of (1.2), the number of glide directions A1 > 4 is even. As has a cyclic ordering, we assume 
gi, ..., gAT to be ordered counter-clockwise, and we set for convenience go := gAr and gAi+i := gi- 
Then, we define the bisectors 


Si ■■= 


gi -k gi+i 


G i = I,..., N. 


\Si + Si+i\ 

We define Ag. as the cone spanned by the bisectors g^_^,g^ surrounding gj, i.e. 

Ag^ := {ag'_i -k ^g- I a, /? > 0} C 

As a consequence, 

Ag = |x G I g is the unique maximizer of g • x for all g G G}- 


( 2 . 1 ) 


( 2 . 2 ) 


The set G induces the crystalline norm || • || given by (1.6). The unit ball of || • || is given by co^. 
The dual norm reads 


X L := maxg • x. 
g&5 


(2.3) 



Figure 3. Example of a set glide directions G satisfying (1.2). The figure also depicts 
several related objects such as: the bisector g 4 , the cone Ag^, and the unit sphere of 
the crystalline norm || • ||. 
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2.3. MMS and EDI. We consider 


defined in (1.8), and some E ^ 


as the state space equipped with the ‘quasi-distance’ D 
n as the energy functional. We note that 


this class of energy functionals includes the energy Ss introduced in Section as a regularization of 
the energy S (1.4) which describes the system of screw dislocations. 


MMS: For a time step r > 0, an initial condition G and a set ^ of glide directions 


satisfying (1.2), we consider the D-MMS given by 


rk+l 


G argmin<h(Z:^,X, r), A: = 0,..., [T/r] — 1, 

xeC"' 


z2 = z°. 


in which the functional is given by 


$ : 


^ 2 )^ X (M 2 )« X ( 0 ,oo) ^ (-00,00], ^>(X,Y,r) := 


D2(X,Y) 
' ^ 


+ E(Y). 


(2.4) 


(2.5) 


The only difference with (1.9) is that we use the more general energy functional E. Existence of 
a solution (Z^)^ to (2.4) is guaranteed by E being bounded from below and lower-semicontinuous. 


Uniqueness does not hold in general. 


For any solution (Z*)^ to (2.4) we define two interpolation curves as mappings from [0, T] to (M"^) 


2\n 


The first one is the piecewise constant interpolant 

7ir{t) := Z^ for t G {{k — l)r, fer]. 


( 2 . 6 ) 


The second one is a De Giorgi interpolant which is defined as follows. For any t G [0, T], let 


:= [t/r] — 1, and note that t G {Tkr,T{kr + 1)]. 

Let 6 := t — rkr G (0, r], and finally 

Z^(f) G argmin<h(Zj:^,X,(i). 
xe( 


(2.7) 

( 2 . 8 ) 

(2.9) 

The main challenge in passing to the limit r —?■ 0 in (|2.6| is that D is not a metric. If it were a 
metric, the techniques in |AGSn81 Ghapter 2 and 3] would apply directly. For this reason, it is useful 
to consider the largest metric d which is smaller than d. It is easy to see that d is induced by the 
crystalline norm || • || on defined in (1.6). We further set 


As (5 < r, the following basic property (see the proof of [AGSn8| Lemma 3.1.2]) holds 

D(Z^,Zr(t)) <D(Z^,Z,(t)) =D(Z^,Z^+i). 


d{x,y) ■■= I|x - y| 


d2(X,Y) := > 


E-'"' 

i=l 


( 2 . 10 ) 


We define the D-MMS, and Z(; analogously to (2.4)-(2.5) by replacing D with D. 


EDI: Now we introduce the energy-dissipation inequality (EDI). To this aim, we first show that (1.3) 
can be written as a (generalised) gradient flow (cf. [M ie H). 


Definition 2.1 (Generalised gradient flow). A triple (Y, E,TZ) is called a generalised gradient flow 
if 

• the state space Y is a smooth manifold; 

• the energy functional E : Y —)• (— 00 , 00 ] is smooth enough for the subdifferential D E(x) G 

Y to be well-defined for all x G Dom E; 

• the dissipation potential 77. : TY —)• [0,oo] is convex, lower semicontinuous, and satisfies 
TZ{x, 0) = 0 for all x G Dom E. 
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The related evolution is given by 

X (x, — D E(x)), in X for a.e. t G (0, T), (2-11) 

where TZ*{x, •) is the Legendre transform of TZ{x, •), and denotes the subdifferential with respect 
to the second argument of TZ*. 

We consider the generalised gradient flow given by the triple iL), where 'L is defined by 

its Legendre transform 

V^WM2^[0,oo), Jlieil*, 




'r^[o,oo), 


( 2 . 12 ) 


2=1 


With this choice, we obtain (1.3) as the evolution (2.11) of the triple Li, T). This is easy to 

see after computing 

= X ... xDV’*(U C 

ro, ifc = 0, 

{(g • ^)g}, if g G ^ the unique maximizer of g • 


DV’*(0 = < 


co{(g'-0g',(g"-6g"}> if {g',g"} = argmaxg-^. 

g&G 


By the basic properties of the Legendre transform, we obtain 

n 


V'(x) = ^||x| 




Z=1 


We remark that by the Euclidean structure of we can identify the tangent space and its dual 

at any X G (M^)” with On the other hand, we do distinct in our notation between particle 

positions X, velocities X, and forces because these objects have a different interpretation. Likewise, 
we use the metric D to measure the distance between two particle configurations, || • || to measure 
velocities, and || • ||* to measure forces. 

The evolution can be written equivalently as a force balance or as a power balance. An overview of 
these descriotions is eiven in |Mie n. By using that E G C'^((M^)"'), |Miel41 Theorem 3.2] provides 
a fourth equivalent description of (2.11), called the EDI. For Li, 4'), the EDI reads 

(2.13) 


E{Z{T)) - E{Z{0)) + / T (Z(t)) + T*( - VE{Z{t))) dt < 0. 


To define a solution concept for curves Z satisfying ( pT^ , we define the space of absolutely 
continuous curves |AGS08I Def. 1.1.1]. For p G [l,oo], we say that X G AC'^’(0, T; (M^)"') if there 
exists an f € LP{0,T) such that 

|X(s) -X(t )|2 < ^ f{r)dr, for all 0 < s < t < T. (2.14) 

By |AGS08[ Rem. 1.1.3], we have for X G ^(7^(0,T; (M^)”) that its derivative is defined almost 
everywhere. 

Definition 2.2 (Solution to EDI). A curve Z : (0,T) 

AC^{0,T; (M2)^) satisfies ( |2T^ . 


(M^)” is a solution to the EDI if Z G 
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The EDI can also be written in terms of the right metric derivative of Z(t) and metric slope 
of E(Zi{t)) (for the precise definition and basic properties, see e.g. |AGS08[ Chapter 1]). For our 
purposes it suffices to define them, with respect to the metric D, respectively as 

sit s-t ^ x->z D(X, Z) 

The right metric derivative and the metric slope with respect to the ‘quasi-distance’ D are defined 
analogously. 

Lemma 12.31 and Lemma 12.41 below contain the standard relation between T-T* and the metric 


(2.16) 


derivatives and slopes. Together they imply that the EDI (2.13) is equivalent to 

r-T 


E{Z{T))-EiZ{0)) + ^ 


'0 


\Z'[%{t) + \dE\UZ{t))dt<0. 


In addition. Lemma 


2.4 


o'- ' ' 

guarantees equality of the slopes with respect to D and D. 


2.4. Basic estimates. In the following two lemmas we prove the relation between the couple T-T* 
and the metric derivatives and slopes when related to D and D. 

Lemma 2.3 (Relation between metric derivative and T). For any Z G AC'^([0,T]; (M^)"'), it holds 

1 , 


-\z'\l{t) = ^[z 


for a.e. t £ (0, T). 


Proof. Since d is induced by the norm || • ||, it holds for a.e. t G (0, T) that 

|2 n 


|zT(*) = i]“E' 

2 = 1 


*(s) -Zi(t)|| 


(s-ty^ 




= 2T Z 


2=1 


□ 


Lemma 2.4 (Equality of slopes). For E G C'^((M^)"'), it holds 

\\dE\l = \\dE\l = ^*{VE). 

Proof. The second equality is given by [AGSOSl Corollary 1.4.5]. For the first equality, we have 
\dE\£) < by the definition of D and D. For the opposite inequality, we observe 

ve;(x)-(x-y) + o(|x-y|2 ) 


|5^Id(X) = limsup 

Y-s>X D(X, Y ) 

Next we are going to construct a particular sequence Y^ —?• X as e 

gi G arg max [g • ViE{X)] , yf := Xi + e [g* • ViE;(X)] g* 
gee 

where VjEl(X) := VxjEl(X) G We note that, by definition of || • ||*, 

gi■V^E{X) = ||ViE(X)||, =:ui. 
Using the explicit sequence (Y^) in ( |2.17 ), we obtain that 

EiLi • ViE;(X) -h o{eui)] 


(2.17) 


0. Let 

i = l,...,n, 


\dE\D(^) > limsup- 
£—^0 




elul,-|-ofeluL) , , , -;-^—rv 

= limsup ' ^ = |u |2 = V2^*(VE;(X)). 

£-s>o e|u|2 


□ 


Now we prove two estimates that are crucial for the proof of Theorem 3.1 
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Lemma 2.5 (Estimate on single time step MMS). Given an energy E such that E G C'^((]R^)"') n 
there exists C > 0 such that for any r > 0, X G we have 

D{X,Y)<Ct, and D{X,Y) < Ct, (2.18) 

where Y,Y G (M^)" are a minimisers of respectively <l>(X,-,r) and ^>(X,-jt). 

Proof. We have 


E (X) > inf 

ze'"'" 


(X,Z) 

2t 


+ E{Z) 


(X,Y) 


2r 


+ EiY) 


and hence, by Lipschitz continuity of E, we obtain (2.18) from 

(X,Y) 


2r 


<E{X)-E (Y) < CD (X, Y). 


Since the proof for D is analogous, we omit it. 

Lemma 2.6 (Bound on directional slope). Let E ^ 
T > 0, andY G argmin<l>(X, -jt). Then 


^)”) bounded from below, X G 


□ 

[)2\n 


D{X,Y) 


> sup 


lim 

Ye-^Y 


E{Y)-E(Y,) 

D{Y,Y,) 


(yf - Xj) // (yi - Xj) for alii = l,...,n 


Remark 2.7. Lemma 2.6 follows directly from the proof of Lemma 3.1.3 from |AGS08) . We interpret 


the right-hand side of the inequality in Lemma 2.6 as a ‘directional slope’. It is clear from the 


dehnition of the slope (2.15) that the directional slope is smaller than or equal to the ‘standard’ slope 
\dE\D{Y). _ 

In fact, |AGS08I Lemma 3.1.3] states that if D were a metric (i.e. satishes the triangle inequality), 
then a stronger inequality than the one in Lemma |2.6| would hold, in which the right-hand side is 
replaced by |9iil|£)(Y). Such a stronger inequality is required in |AGS08| to prove convergence of the 
MMS to the EDI. 

Since D is not a metric, we cannot apply this stronger inequality. It is not hard to construct an 
example in which \dE\jo(Y) is indeed larger than the ‘directional’ slope. The bulk of our proof of 
the main result (Theorem [3.1 [) concerns an alternative argument to |AGS08I Lemma 3.1.3]. 


3. Convergence of the D-MMS to the EDI 
We now prove the main result of the paper. 

Theorem 3.1. Let T > 0 be an end time, G (M^)” an initial condition, and E G (^^((M^)"^) be 
an energy such that E and YE are bounded and uniformly continuous. For any time step t > 0, let 
(Z:^)fcgN be a solution to the D-MMS (2.4) and Z,- : [0,T] —)• (M^)” the corresponding step function 
defined in (2.6). Then, along a subsequence o/r —>■ 0, the curves Zr{t) converge pointwise for all 


t G [0,T] to an absolutely continuous curve Z : [0, T] 


p2\n 


)”■ which satisfies the EDI (2.16). 


Proof. We split the proof in four steps. In Step 1 we use a rehned version of Ascoli-Arzela to identify 
the limiting curve Z G AC‘^{[0,T]] (M^)”) for a subsequence of Z^ as r —>• 0. Step 2 is a consequence 


of Step 1, in which we precisely state in what sense the interpolants Z^ and Z^ (2.8) converge 


Step 3 is the main novelty of this paper, in which we provide an alternative argument to |AGS08l 
Lemma 3.1.3] which allows us in Step 4 to pass to the limit as r —?■ 0 to obtain that term in the 
EDI which is related to the slope \dE\j^. In Remark 2.7 we discuss the relation between |AGS08l 
Lemma 3.1.3] and the weaker result in Lemma 2.6 
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Step 1: Compactness of the step functions Z,-. 

We prove that Z,- converges point-wise in time along a subsequence to some 

|AGSn8[ Proposition 3.3.1] guarantees the existence of Z G C([0, Tj; provided that there exists 

a symmetric function ui G [0,T]^ —>■ [0, oo) satisfying 

limsup (Zt-(s), ZT-(t)) 0. (3.1) 

T^O 


We prove (3.1) by estimating D (Zt-(s), ZT-(t)) for arbitrary G [0,T], where s <t without loss of 
generality. Let K = \s/t\ and L = [f/'r], and note that s G {{K — l)r, Kt\, t G {{L — 1)t, Lrj. Then 
by the definition of Z^ (2.6), we have Zt-(s) = and ZT-(t) = Z:(i, and we estimate 

L-l 

rk+l 


b (Zr{s), Zr{t)) = b (zf, z^)<J2d (z';, z 


k=K 


We continue the estimate by using Lemma 2.5 to obtain 

L-l 


b {Zr{s),Zr{t)) < ^ Ct = C{L - K)t < CT{\t/T] - [s/t]) 


k=K 

< Ct ^ + 1 ] = C{t - s) + Ct, 


(3.2) 


from which we conclude that (3.1) is satisfied for Lo{s,t) = C\t — s|. In the sequel, we proceed with 
the subsequence provided by |AGS08I Proposition 3.3.1] without changing notation. 

We prove that Z is absolutely continuous by showing that (2.14) holds for some constant function 
/. From ZT-(t) —)■ Z{t) for a.e. t G (0,T) and (|3.2|), we deduce that 


b ^Z(s), Z{t)J = lirn L) (Zr{s), Zr(t)'j < J C, for a.e. 0 < s < t < T. (3.3) 

Since Z is continuous, we conclude that (3.3) holds for all 0 < s < t < T. 

Step 2: Convergence of interpolants Z^ and Z^ (2.8). 

We fix t G [0,T], and let kr = It/T] — 1 as in (2.7). The compactness resu lt in Step 1 implies directly 
that Z1(Z^^_|_]^, Z(t)) —)• 0 pointwise in t as r —)• 0. Then from Lemma 2.5 and (2.9) we also have 

b{z^^-,zit))+b{zlit),zit)) ^0. 

As a result of this and \/E being continuous, we have that VFl evaluated at Z^^, and Z(-(t) 

converge to VE{Z{t)) as r —)• 0. 

Step 3: Convergence of directional derivative to slope. 

We fix an arbitrary t G [0,r] and set kr as in (2.7). Let 

o’" = {g[,... ,g^} G be such that gf//[Z(((t) - Z^-]. for alH = 1,... ,n. (3.4) 

We prove in this step that 


[g[ • V.EiZlit))]' ^ \dE\l{Z{t)). (3.5) 

i=l 
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From the proof of Lemma 2.4 and the definition of 4^* ( |2.12[ ), it follows that (3.5) is implied by the 
claim 


[gf • ViF;(Z^(t))] max [g • ViF;(Z(t))] , for all z = 1,..., n. (3.6) 

g65 


For proving (3.6), we set a* := and g^ as a minimiser of —g • a*. If a* = 0, then (3.6) 

follows directly from Step 2 . 

The main part of the proof of (3.6) for aj 7 ^ 0 is to characterise G^. We recall from (2.8) that 


Z^(t) G argmin 


y^D\z';^,X) + E{X) 


(3.7) 


xe(R2)’^ 

where 5 = t — rkr G (0,r]. Since D is only finite on MQ, we can restrict the minim isati on over 
X = (xi,... ,x„)^ to Xj = + ctjgi with cii > 0 and gj G Q. Thanks to Lemma 


2.5 


we can 


assume that ai < C6. At the same time, we expand E{'K) in a Taylor series around Z^^. By these 
arguments, we rewrite the minimisation problem in (3.7) as 


E{Z^'^) + min min 

giSS OLi>0 

1=1 


4«? + aigi-V,.E(Z^) + o(,5) 

AO 


(3.8) 


We recall from (3.7) that gf is a minimiser of this minimisation problem. We characterise this 
minimiser by solving (3.7) first under the assumption that o{S) = 0. Minimising (3.8) over a, > 0 
yields 

7^(Z^) - - y; (min g, • V.i^(Z^)) . (3.9) 

i=\ 

We treat the minimisation within parentheses for each i = 1,... ,n separately. For convenience we 
drop the index i in our notation whenever possible. We split characterizing g^ in two cases: (i) —g • a 
has a unique minimiser g, and (ii) —g • a has exactly two minimiser g^ and g^. 

In case (i), it holds that a is an element of the cone Ag (2.1). Hence, there exists an r > 0 such 
that B(a,r) CC Ag. Then, by Step 2, it holds that ViE{Z^^) G B{a,r) for all r small enough. 
Hence, g is the unique minimiser of g • Vi£'(Z^^), and, furthermore, there exist a c > 0 independent 
of r such that 

{ min (g - g) • ViFl(Z^^)) > c > 0. 

Hence, the minimiser of g • Vjiil(Z^^) is stable under perturbations of the form o(l), and thus gj 
minimises the term in square brackets in (3.8). Hence, gf = gj, from which we conclude that (3.6) 
holds in case (i). 

In case (ii), it holds that a is an element of the cone A spanned by g^ and g^. An analogous 
argument as used in case (i) implies that g • Vjiil(Z^^) is minimised by g^ or g^, and, furthermore, 
there exist a c > 0 independent of r such that 


( 


min (g — g'^) • ViE{Z ^'^)) > c > 0 , for j = 1 , 2 . 
g6d\{gLg2} 


Hence, either g^ or g"^ minimises the term in square brackets in (3.7). This implies g^ G {g^,g^} 
and (3.6) follows from the observation that 

r •v,i?(zr(t)) 




> g-' • a = mm g • a, 
g65 


for j = 1 , 2 . 


Step 4: Z satisfies the EDI. 


To simplify the proof, we assume that Z^ and Z\. are uniquely defined by (2.4) and (2.8). In the 
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general case we should take into account the supremum and the infimum of over all 

possible instances of Z^(t). This is done in |AGSn8| Chapter 3], in which all equalities below become 
inequalities. 

Following the lines of |AGSn8| Theorem 3.1.4], we obtain for a fixed time step r > 0 that 


D 2 ( 2 ^, z^+i) rAfe+i) Z^) 


2r 


+ 


Summation over k results in 

[T/rl-l 




2{t — rkY 

/■rr/r]T 


dt = F;(Z^) - E{Z^+^ 


for A: = 0,..., [T/r] — 1. 


fc=o 


2t 


2{t- [t/rj rf 


dt = E{Z^)-E{z']^/^^). (3.10) 


Next we show how to pass to the limit in the terms of (3.10) as r —)• 0. We set kr{t) = [f/r] — 1 
as in (2.7). By Step 2, the right-hand side of (3.10) converges to the difference between the energy 
evaluated at Z(T) and Z^. Regarding the first term in the left-hand side of (3.10), we use the estimate 
|AGS08l (3.3.10)] to obtain that 

r^2(rjk ’7k+l\ 


r) 2 ( 7 k rj 

liminfr V ^ 

^ ‘2'r2 


k=0 


2r2 


= liminf 
r^O 


qr/rlr £)2 (zGW 1 


2r2 


1 f'^ 


dt> - 

Z' 

~ 2 Jo 



D 


(t) dt. 


We prove the convergence of the second term in the left-hand side of (3.10) as follows. We fi x t G [0, T] 
and r > 0. We set as in 


and Ui := gl • ViE{Ztj.{t)). We obtain from Lemma 2.6 that 


---—-- > hmsup^ Veigf ■ViE{Z^{t)) = hmsup 

t — Tkr[t) e_>0 ^ e-s>0 

|2 


e ■ u 


^2 


= |U|2- 


In Step 3 we prove that |u |2 —)• |clFl|^(Z(f)) as r —)• 0. Hence, by applying Fatou’s Lemma, we obtain 


lim inf - 

T —>-0 2 


1 D\Zl{t),z\^'^^) 1 

(t - [t/rj r)2 - 2 , 


\dE\l{Z{t))dt. 


Combining the results above, we obtain after passing to the limit r —)• 0 in (|3.10) 


i'/|2 


1 


'0 


' 2 

from which we conclude by Lemma 


{t)dt+- / \dE\%{m)dt<E{Z^)-E{Z{T)), 


'0 


2.3 


and Lemma 2.4 that Z satisfies the EDI. 


□ 


4. Generalisations of Theorem 13.11 


Here we discuss in more detail two of the three remarks on Theorem 13.11 that are mentioned in 
Section 11.51 


4.1. Extension to higher dimensions. Theorem 3.1 can be generalised to higher dimensions for 
the particle positions. We show this by considering n particles with positions Zj G with d > 2. 
Since the setting is analogous to the two-dimensional scenario, we keep the same notation, and only 
mention the important differences. 

The main difference in higher dimensions is that we consider the set of glide directions given by 

g := {gl,... ,gAr} C 

which satisfies the two properties 

gG^ -g£g, and span^ = M‘^. 
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The related crystalline norm reads 


N 


k=l 


N 


x|| := min < ak Ofc > 0 such that E ^ I ■ 


k=l 


Next we show that the dual norm can be characterised by 

X • y 


X L, := max 


yeR'^MO} ||y| 


= max g • X. 

g66 


(4.1) 


(4.2) 


Since the maximum is taken over a smaller space in the right-hand side of (4.2), it suffices to prove 


that it is larger than or equal to the left-hand side. To this aim, let y be a maximiser of the left-hand 


side, and let > 0 as in (4.1). We conclude by estimating 

N X -1 iV 


XL = 


X- y 


k=l 


y^Okgk ■ X < maxg • X. 


In terms of the crystalline norm and its dual we define the distances D and D and the related T and 
T* analogously to the two-dimensional setting. With these objects, the MMS (2.4) and EDI (2.13) 
are defined analogously. 


It is readily checked that the basic estimates in Section 2.4 and most steps in the proof of Theorem 


3.1 hold in the d-dimensional case by analogous arguments. The only part where the extension to d 
requires a modification is for the characterisation of the minimiser of the term within parenthesis in 


(3.9) in Step 3 in the proof of Theorem 3.1 The remainder of this section describes this modification. 


First, we define the cone Ag as an extension of (2.1) by 


Ag := 


X e 


g • X > max g • X 

g6d\{g} 


(4.3) 


We remark that {Ag n ^}gee describes the Voronoi tessellation of S"^ ^ with respect to Q. 


d-l 


With (4.3) it is easy to see that the characterisation of the minimiser in (3.9) can be done with an 


analogous argument when a G Ag for some g G If a ^ UgggAg, then 

argmaxg • a > 2, 
g65 

while in the two-dimensional scenario the left-hand side equals 2. In any case, the argument in Step 3 
also holds when more than 2 glide directions are considered. 


4.2. Annihilation of dislocations in a variational framework. With Theorem 3.1 we proved 
that the time-continuous model in |CG99j can be obtained as the limit of the time-discrete D-MMS 
schemes when the time step converges to 0. We would like to address an open problem regarding 
annihilation of dislocations. 


By regularizing the energy to prevent the minimising-movement scheme in (1.9) to jump in the 


first time step to a state in which the energy equals —oo, we remove the possibility for dislocations to 
annihilate or to be absorbed at dD. This regularization of the energy is therefore artificial. Indeed, the 


evolution defined by (1.3) can easily be modified to allow for annihilation; whenever two dislocations 


annihilate (or one dislocation leaves the domain D), the evolution can be restarted by removing the 
annihilated dislocations from the equation, and taking the current positions of the other dislocations 
as the new initial condition. This raises the following question: is it possible to modify the MMS to 
describe such dynamics which allow for annihilation in a variational framework? 
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In |ADLGP14| the minimum in (1.9) is taken over a r-independent neighbourhood M of With 
this strategy, the energy does not need to be regularized, but the MMS breaks down whenever the 
difference between any pair of opposite dislocations is in N. 

One idea to make this approach work, is to make A/" dependent on r, such that it shrinks to a 
singleton as r —?• 0, but at a slower rate than r (e.g., \/r). 

The main challenge is to treat the time step in which annihilation takes place. From the procedure 
above, it seems reasonable to restart the MMS with a new energy (consisting of fewer dislocations) 
while keeping fixed the positions of the dislocations which were not annihilated. We expect that 
passing to the limit r —?• 0, if possible, requires serious modifications to the argument in [AGSOSI 
Ghapter 2 and 3]. 
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